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We present a quantum mechanical description of parametric down- conversion and phase-matching 
of Bloch-waves in non-linear photonic crystals. We discuss the theory in one-dimensional Bragg 
structures giving a recipe for calculating the down-converted emission strength and direction. We 
exemplify the discussion by making explicit analytical predictions for the emission amplitude and 
direction from a one-dimensional structure that consists of alternating layers of Alo.4Gao.6As and 
Air. We show that the emission is suitable for the extraction of polarization-entangled photons. 
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I. INTRODUCTION 

Entangled photon pairs play a central role in both fun- 
damental tests of quantum mechanics and in the imple- 
mentation of quantum information theory protocols 
They are an appealing resource for quantum communica- 
tion since they propagate easily over long distances with 
relatively little interaction with the environment. Fur- 
thermore they have been proposed as a resource for all- 
optical quantum computation Q. 

A popular method to produce entangled photons is 
by parametric down-conversion in naturally birefringent 
non-linear crystals. In this process 'pump' photons enter- 
ing a crystal decay or 'down- convert' into photon pairs. 
The role of the non-linearity is to mediate the interaction 
whereas the role of the birefringence is to ensure that the 
process is phase-matched, i.e. that the amplitudes for the 
down-conversion process at different points in the crystal 
constructively interfere. If a particular emission geome- 
try can be achieved, photons emerging in a specific pair 
of directions will be entangled in polarization . 

In a recent paper, the authors proposed the use of 
non-linear photonic crystals as a source of polarization- 
entangled photons 3- The scheme harnesses the 
higher nonlinearity of semiconductor materials (e.g. 

XgL = 200pm/V i c.f. xtlo = 2.2pm/V 0) to 
mediate the down-conversion and proposes to use the 
photonic crystal geometry to phase-match the emission. 
The scheme has the potential of both increasing the effi- 
ciency of the process and providing an entangled-photon 
source that is more amenable to integration on opti- 
cal chips. The scheme is fundamentally different from 
schemes for quasi phase-matching in periodically-poled 
materials, where only the is modulated periodically. 

One-dimensional photonic crystals had been consid- 
ered before for classical frequency conversion in the limit 
that the optical wavelength is much larger than the pe- 
riodicity [3, |g, lij • Two-dimensional non-linear photonic 
crystals have also been considered for classical frequency 
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conversion 0. Following the appearance of our pro- 
posal, Reference El used semi-classical coupled mode 
theory to calculate co-linear frequency down-conversion 
efficiency in one-dimensional structures using numerical 
calculations. 

In this article we present the theory that underlies our 
proposal. We use the Bloch-wave formalism to discuss 
the quantum mechanical down-conversion process in non- 
linear photonic crystals and show how to calculate the 
strength and direction of the down-converted emission. 
The theory works for photonic crystals of all dimensions. 
We apply the theory to one-dimensional Bragg struc- 
tures, performing calculations of the down-conversion 
emission analytically. By plotting the calculated emis- 
sion from a structure that consists of alternating layers of 
Alo.4Gao.6As and Air, we show explicitly that entangled 
photon pairs can be generated in a realistic structure. 

The present article is structured as follows: sections ITTI 
and IIIII discuss Bloch waves and their quantization in 
linear photonic crystals following essentially the work of 
Caticha and Caticha |l2|. In section HVl we derive an ex- 
pression for the interaction Hamiltonian and the phase- 
matching function that govern the down-conversion pro- 
cess in non-linear photonic crystals. In section^we sum- 
marize the discussion of Bloch waves in Bragg structures 
by Yariv and Yeh 13] extending it to obtain expressions 
for the Bloch-wave Fourier components. Section con- 
sists of a detailed discussion of the phase-matching prob- 
lem in Bragg structures, giving a recipe for establishing 
the emission amplitude and direction from a given struc- 
ture. The discussion centres on the application of our 
method to an example Alo.4Gao.6As/Air structure. Fi- 
nally we conclude and discuss extensions of the present 
work in section IVTTl 



II. WAVE PROPAGATION INSIDE A LINEAR 
PHOTONIC CRYSTAL 

A photonic crystal is a material with a periodic vari- 
ation in the index of refraction or dielectric constant. 
A periodic dielectric constant e(r) = e(r + A), is seen 
as a periodic potential by the electric field. The eigen- 
solutions of Maxwell's equations in a medium with pe- 
riodic dielectric must therefore take the form of Bloch 
waves. We can thus write the following expression for 
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the four- vector potential inside the medium: 
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K represents the Bloch momentum (and will be taken 
to lie within the first Brillouin zone), the index A runs 
over the two polarizations and the various branches of 
the dispersion relation, G represents a reciprocal lattice 
vector and \i runs from to 3. The bold face represents a 
Cartesian vector and only the real part of is of physi- 
cal significance. The sum defines a periodic envelope-field 
a K \( r ) = a K A( r +^) W ^ n Fourier coefficients qk,a(G). 
There are two natural generalizations of the Coulomb 
gauge: e(r)V • A = and V • (e(r)A) = that reduce 
to the uniform dielectric Coulomb gauge V • A = 0. In 
the gauge e(r)V • A = 0, Maxwell's equations are given 
by 
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Solving these equations with the ansatz (Eq. [Q amounts 
to finding the dispersion relation between uo and K 
and an expression for A (r) or its Fourier coefficients 
d^(G). The most striking consequence of the pres- 
ence of a periodic potential is the formation of frequency 
regions, known as 'stop-bands' in which no propagating 
solution exists. Close to these regions, the dispersion is 
strongly modified and c^k,a (j*) takes the form of a stand- 
ing wave. Away from these regions «k,aW recovers its 
plane-wave form but the dispersion relation can differ 
considerably from that in a uniform dielectric in a way 
that can be tuned by changing the geometry and strength 
of the periodic modulation. The tunable dispersion rela- 
tion will play a crucial role in this paper. 

The electric and magnetic field can be derived from 
the vector potential in the usual way: 

E k ,a = -V0k,a - <9 £ A K ,A 
B k ,a = V x A k ,a 

Where k ,a = 4,a and A = (^k,a^k,a^k,a)> g^- 
ing: 

E k ,aM) = e-^ K -™'^^s K ,A(G)^ G - r (4) 

G 

B k ,aM) = e-^ K -™'^^/3 K ,A(G)e^- r (5) 

G 

where the Fourier coefficients sk,a(G) and /3k, a (G) are 
given by: 

sk,a(G) = -i(K-G)^ K?A (G)-io;aK,A(G) (6) 
/3 k ,a(G) = iGx d K ,A(G) (7) 

The symbol for the Fourier coefficients of the electric field 
e should not be confused with the dielectric constant e(r). 



III. FIELD QUANTIZATION INSIDE A LINEAR 
PHOTONIC CRYSTAL 



Unlike sum frequency generation, parametric down- 
conversion does not occur classically and thus is a truly 
quantum-mechanical phenomenon. It is thus most natu- 
ral to discuss the problem in the language of the quan- 
tized electromagnetic field. The procedure for the quan- 
tization of the electro-magnetic field in a medium with 
non-uniform dielectric differs from that in a uniform di- 
electric in that the equations of motion for the poten- 
tial (Eqns.[3E|) involve derivatives of the dielectric func- 
tion. The two natural generalizations of the Coulomb 
gauge: e(r)V • A = and V • (e(r)A) = lead to dif- 
ferent quantization procedures. Here we summarize the 
results of Caticha and Caticha who used the gauge 
e(r)V • A = to quantize the electro-magnetic field in a 
medium with periodic dielectric. They showed that the 
Hamiltonian is diagonal in the Bloch-wave basis and the 
creation operators for the field satisfy suitably modified 
commutation relations. The quantized field operator is 
given by: 



A 



d 3 K 
(2tt) 3 



a(K,X)A^ x (r,t)+h.c. (8) 



Where A (r, i) is given by Eq.[l] The Hamiltonian can 
be expressed in the usual form: 



H 



/H 3 K 
^3^ K ,Aat(K,A)a(K,A) 



and the creation and annihilation operators a and a) sat- 
isfy the following commutation relations: 

[o(K,A),ot(K',A')] = (2tt) 3 ,5( 3 )(K-K')<5a,a' (9) 
[a(K,A),a(K',A')] = (10) 
[ot(K,A),ot(K',A')] = (11) 

The electric and magnetic field operators derived from 
Eq. [S| are given by: 



E(r,i) 
B(r,t) 



E 



d 3 K 
(2tt) 3 

d 3 K 
(2^) 3 



a(K,A)E K ,A(r,£)+/i.c. (12) 



a(K, A)Bk,a(f, t) + h.c. (13) 



with EK,A( r ?£) an d BK^fr^) given by Eqns. 21 and 
We will now use these results to derive an expression 
for the quantum interaction Hamiltonian and the phase- 
matching function. 



IV. NON-LINEAR PHOTONIC CRYSTAL 
INTERACTION HAMILTONIAN AND 
PHASE-MATCHING 

To derive the interaction Hamiltonian for the 

quantized electromagnetic field, we proceed as in the case 
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of non-linear optical crystals with uniform dielectric and 



(2) 



Starting with the expression for the classical 



X 

interaction Hamiltonian: 



/dvi, 



(where is the second order susceptibility tensor), in- 
serting the expression for the quantized electric Bloch- 
field E(r, 0) (Eq. IT2|) . making the rotating wave approx- 
imation and labelling the three interacting Bloch modes 
by p, 1,2, we obtain the following quantum interaction 
Hamiltonian: 



p,l,2 



dV 



^oxSlW (r)Ei(r)E^r)a p a\al + h.c.(14) 



tAt 



Where ^,1,2 is short for £k p ,i, 2 ,a p ,i, 2 ( r > 0) and S,i,2 is 
short for a(K p> i 5 2, A^i^)- The interaction can be seen 
to mediate two basic processes: one in which a p photon 
down-converts into photons 1 and 2 and (h.c.) in which 
photons 1 and 2 up-convert to photon p. 



To calculate the time evolution of the field under this 
interaction Hamiltonian, we switch to the interaction pic- 
ture, and evaluate the first term in the Dyson series ex- 
pansion for the time evolution of an initial state |^o)- 
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For given states p, 1,2, the amplitude for parametric 
down-conversion into modes 1 and 2 is proportional to 
the phase-matching function <3>(p, 1,2), given by: 



$(p, 1, 2) = e I dV X %lE?(r)Ei(r)E*(r) (15) 



It is instructive to substitute the Fourier expansion 
of E p (r), ^i(r), E2(r) (Eq. 0} and the Fourier expansion 

(2) (2) 

of Xijk( r ) ( w ith Fourier coefficients denoted by X^/~(G)) 
into the phase-matching function (Eq. IT5|) to obtain: 



*(p,l,2) 



G x ,G P ,Gi ,G: 



eoxgi(G x )e; i (G P ) £ l(G 1 ) £ 1(G 2 ) ( 5( 3 )(K p 



Ki — Ko 



Gi + G2 — G^ 



(16) 



which makes the conservation of Bloch quasi-momentum 
manifest. 

Equations El and El are the main result of this sec- 
tion. To first order in a pump photon will down- 
convert into a superposition of all Bloch-wave pairs that 
satisfy conservation of energy (c^i+c^ = w p ) and of Bloch 
quasi-momentum, with an amplitude <£(p, 1,2) given by 
Eqn. El or Eqn. El An important difference with the 
phenomenon of up-conversion is the fact that down- 
conversion involves all modes that phase-match, whereas 
up-conversion is more constrained: given two photons 
to up-convert, there is typically only one mode they 
can up-convert to. A parallel between phase-matching 
in non-linear crystals with uniform dielectric and phase- 
matching in photonic crystals is drawn in Table [I] which 
contains a dictionary for the main concepts. 

The vector equation that expresses the conservation of 



Bloch quasi-momentum: 

Kx + K 2 = K p + G x + Gi + G 2 - G p (17) 

shall be referred to as the phase-matching equation. For a 
given K p , solving the equation corresponds to finding the 
intersection between the dispersion surfaces |Ki(fci,o;i)| 
and |K2(^2,^2)| for photons 1 and 2, centred on the ori- 
gin and on \JZ p (k p ,(Ji + ^2) | where k represents a unit 
vector that points in the same direction as K. To solve 
the equation it is therefore necessary to compute the dis- 
persion surfaces of photons with frequencies ui,U2 and 

For photons that satisfy the phase-matching equation, 
the amplitude of the process is proportional to the over- 
lap of the waves in the non-linear medium. As we shall 
see in later sections, this is most easily calculated us- 
ing Eq. El i n cases where the light field has standing 
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Natural non-linear crystals 


Non-linear photonic crystals 


Plane waves 


Bloch waves 


E A (r) = e* kr 


E K x(r) = e* Kr V e K A(G)e^" r 
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Momentum conservation 


Quasi-momentum conservation 


ki + k 2 = k p 


Kl + K 2 = Kp + (or 
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birefringence 


form birefringence 


k(u;, k) = n((jj, k)/c 


k(a;,k) = K(cj,k) 


Tensor properties x^ 


Tensor properties x^ an d 


determine amplitude 


Fourier coefficients 




determine amplitude 



TABLE I: Phase-matching dictionary that shows how the 
main ideas of phase-matching in non-linear optical crystals 
translate to non-linear photonic crystals. 



wave character and using Eq.[HDwhen it has the charac- 
ter of a propagating wave. It is therefore useful to keep 
both these expressions in mind when solving the phase- 
matching problem in a given structure. 

In deriving Eq. [03 the photonic crystal (interaction) 
volume was taken to be infinitely large. This will of 
course not be the case for real photonic crystals. For 
finite-dimensional crystals the delta functions embodying 
the conservation of momentum become sine functions: 

sin [(K p - Kx - K 2 + G x + d + G 2 - G^U] 



(K p — Ki — K 2 + G x + Gi 



where Li is the length of the photonic crystal in the ith di- 
rection. This corresponds to a certain amount of allowed 
momentum mismatch and, as will be seen in section IVll 
can have significant consequences. 

In one and two dimensional photonic crystals, there 
are two main mechanisms by which the dispersion sur- 
faces are modified by the presence of the crystal: form 
birefringence 7] and geometric dispersion. Although the 
two effects are not entirely independent of each other, 
their origin is physically distinct. Form birefringence is 
the difference in the dispersion surfaces of Bloch waves 
that have different polarizations and arises from the dif- 
ferent boundary conditions at the interfaces in the pho- 
tonic crystal. In the long wavelength limit (A ^> A), this 
type of dispersion is the dominating one. The geomet- 
ric dispersion is induced by the presence of the periodic 
potential and appears as the only type of dispersion for 
waves propagating in directions for which the boundary 
conditions do not break the symmetry between polar- 
izations. In three dimensional structures the problem is 
more complicated since the de-coupling between direc- 
tion and polarization does not readily occur. 

The derivations so far apply in all dimensions, how- 
ever from here on we will restrict our attention to one 
dimensional (Bragg) structures and discuss the analyti- 
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FIG. 1: Illustration of a one dimensional photonic crystal 
composed of alternating layers of materials a and b of thick- 
ness a and b having refractive indices n a and n&. The struc- 
ture is periodic with period A = a + b. The axis of symmetry 
is taken to be the z axis. There are two types of propagat- 
ing polarization eigen-modes: Transverse Electric (TE) and 
Transverse Magnetic (TM). TE(TM) waves have the elec- 
tric (magnetic) field vector lying in a plane parallel to the 
interfaces between the materials. 



cal solution of the down-conversion phase-matching prob- 
lem in detail with a view to generating entangled photon 
pairs. The problem in the two dimensional case is sim- 
ilar, however calculation of dispersion relations in two- 
dimensional structures has so far only been approached 
numerically. For a discussion of phase-matching in two- 
dimensional photonic crystals see refs 0. Ref [D| 
discusses the problem in a structure having a two- 
dimensional periodic x^ 2 \ but uniform linear dielectric. 



BLOCH WAVES IN (ONE DIMENSIONAL) 
BRAGG STRUCTURES 



In this section we summarize the results of Yariv and 
Yeh |l3|, describing the Bloch- wave solutions that prop- 
agate in one-dimensional photonic crystals (Bragg struc- 
tures). We extend their results by obtaining expressions 
for the Bloch-wave Fourier coefficients and discuss the 
form of the Bloch-waves in detail, in order to gain intu- 
ition that will be needed in the discussion of section IVTl 

Figure [I] shows the basic structure of a one dimensional 
photonic crystal consisting of a series of alternating layers 
of materials a and b. We adopt Cartesian coordinates 
with the z axis aligned with the axis of symmetry. Inside 
each layer the field can be expressed as a superposition 
of forward and backward propagating plane waves, with 
wave vectors k having magnitude |k| = ^ (where n 
is the refractive index in the layer) and component k^ b 
along the z axis. Translational symmetry in the plane 
perpendicular to the z axis implies, through conservation 
of momentum, that the components of the wave- vectors 
parallel to the planes (k x and k y ) are equal across the 
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boundaries, thus the z components are given by: 



(18) 



where k 2 



k% + ky and n a ^ are the refractive indices 



of materials a and b. For waves traveling at an angle to 
the z direction, there is a natural choice for the polar- 
ization basis: Transverse Electric (TE) with the electric 
field pointing out of the plane defined by the wave- vector 
and the z axis and Transverse Magnetic (TM), with the 
electric field lying in the plane (see Figure 0). 

We can thus write down the following expressions for 
the electric field in the n th layer: 



Ea A(ujt—k v i 
n,TE(M) = e 



( «n) 



TE(M) e 



-ik z (z — nA) 



(a n £^)te(m) 



Jk a z (z-nA) 



)(19) 



E^TE(M)=^ t ""^ ) ( (^)TE(M)e-^ ( ^ A) + 

(^^-)TE(M)e^ (2 - nA) ).(20) 

where (£± 6 )te(m) are polarization unit vectors defined in 
Figure E] and we have chosen y as the off- axis direction 
for convenience. Matching the fields at the boundaries we 
obtain relations between the amplitudes in neighbouring 
slabs, which can be further reduced to a relation between 
coefficients in the same material in neighbouring cells: 



A B 
C D 



The coefficients for TE and TM waves are 0: 

i tk\ 
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v z ' 2\nlk\ n\k a z ) v z . 



ik a n 



i /nlkz 
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n\k% 



) sin(fe^) 



With C TE /TM = ^TE/TM' ^TE/TM = ^TE/TM* ^he 

eigen-modes of propagation can then be obtained by find- 
ing the eigen- vectors and eigen- values of the transfer ma- 
trix. The right-moving eigen- value is e lKzA with: 



K™/™(fc ||) c) = ^cos- 1 



^TE/TM + ^TE/TM 



(21) 



where K z is the z component of K = (K x = k x ,K y = 
ky,K z ) and the corresponding eigen- vector is: 



a TE/TM = e -inK l z 



b TE/T 



TE/TM 



J-TE/TM 



(22) 
123) 



the coefficients are related to the a^s via (b+,b n ) 
M(a+,a~) where M is given by: 
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(hl±]lD p ia(k^-k b z ) (k h z -k a z ) ia ( k *+ k b) 
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2kl e ^ b e 



2k' 



(24) 



(nlkl+njkl) r ig(kt-k1) (n 2 b k a z -n 2 a kl) i a <k a _+k1) 

1n a n b k\ c 2n a n b k\ c 

jnlK-nlk*) ia(k-+k b .) (njk h ,+nlK) ia ( k "-k b ) 



(25) 

To relate these expressions to the Bloch-wave expres- 
sions in sections |H] and IIIII we need to evaluate the 
Fourier transform of the expressions for the electric 
field (Eqns. EH and 120)1 and compare it to Eq. ^ (with 
G = njf-z). This is done in Appendix 1X1 where we de- 
rive the following expression for sk,a : 

2tt. 



a ^ sm[(K z - K - n% )f] ^_ l{Kz _ K _ n ^^ a 



V2^ ° (K z -k*-n%)* 
a „_sm[(K z + K - n^ e _ i(Ka+k a_ n3jt)iia 



V2^ ° {K z + k«-n^)% 



A A»*T 2 y^ 



(K z -k*-n%)% 

b L _ Sin[(^ z + k b z - »X)|] )(a+4)|6 



V2^° {K z + k» z -n^)\ 

The expressions for the dispersion relations 
|Kp 5 i 5 2(i^p,i,2,^p,i,2)|, the field amplitude E p? i 52 (r), 

and its Fourier coefficients £K p ,i, 2 ,Ap,i, 2 ( n x) can ^ e 
substituted into Equations ^| and EH to give analytical 
expressions for the phase-matched emission from any 
given (one dimensional) non-linear Bragg structure. The 
expressions are however not immediately intuitive. We 
therefore proceed to consider the results of this section 
in various situations to gain some intuition on how the 
Bloch-wave properties depend on the geometry of the 
crystal and on the properties of the materials. 

Central to an intuitive understanding of Bloch- waves in 
a given structure is the band diagram for that structure. 
Figure[2|shows a band diagram for a structure of alternat- 
ing layers of material a with refractive index n a = 1 and 
material b with refractive index n& = 5 (the values were 
chosen for illustrative purposes). Note that the labels a, b 
shall also be used to indicate the thicknesses of materials 
a, b. The fill fraction a/ A of material a is \. For sim- 
plicity, natural material dispersion was neglected; it will 
be included in the section that follows. The gray (white) 
areas correspond to combinations of uo and k\\ that prop- 
agate(do not propagate). A given combination of uo and 
fc|| will propagate if the corresponding K T e/tm(&|| > °°) is 



6 




1 0.5 0.5 1 

Parallel wavevector k,, Or/A) 

FIG. 2: Band diagram for a structure of alternating layers of 
material a with refractive index n a = 1 and material b with 
rib = 5. The fill fraction of material a: a/ A is ~. For simplic- 
ity, natural dispersion and absorption were neglected. Under 
this simplification, the diagram can be drawn for frequencies 
and parallel components of k expressed as multiples of ^ and 
J respectively. The points labelled UJ1...5 correspond to the 
Bloch-wave plots of Figure El and dispersion surface plots of 
Figure 0] 



real. This will be the case if \\(A + D)\ < 1. For TE 
waves \{A + D) is given by: 



cos(k z a + k b z b) 



1 



1 /J4 kl_\ 



whereas for TM waves it is given by: 

1^ k}i 



cos(k a z a + k b z b)- 



n; 



! h,b 



sm(k z a) sin(k b z b) 



sin(/c» sm(k b z b). 



It can be seen that the extent to which K Z A deviates from 
the simple linear relation K Z A = k z a + k b z b depends in a 
non-trivial way on the ratios of the refractive indices and 
on the width of the layers compared to the wavelength 
inside the medium. The term in the square brackets is 
independent of the width of the layers and is a measure of 
the strength of the modification of the dispersion relation 
(and thus the width of the bands). The trigonometric 
terms in turn depend on the layer widths in a simple 
way and have the primary function of determining the 
positions of the bands. It is instructive to consider the 
band-structure for fc|| = of a crystal of fixed optical 
periodicity / = n a a + n^b and varying optical fill fraction 
/ = n a a/l. Measuring the free space wavelength as a 
fraction of I: Af ree S p. = x x I we obtain the following 
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FIG. 3: (Left:) Field strength and (Right:) Fourier compo- 
nents of degenerate TE/TM Bloch waves propagating along 
the axis of symmetry of the illustrative structure with band 
diagram shown in Figure The wave parameters can be read 
off Figure [21 Close to the band edge (002,(^3) the counter- 
propagating components grow to make the wave increasingly 
like a standing wave with the field concentrated in the mate- 
rial of lower or higher refractive index depending on whether 
the frequency is close to the bottom or the top of a band- gap. 
As the wavelength gets shorter, the waves deviate more and 
more from a plane wave even close to the centre of a band 
(w 5 ). 



expression for (A + D): 



/2tt\ 
COS l-J + 



2\n b n a 



. /2tt/\ . /2tt(1- 
sm sin 

V x J V x 



/) 



It can be seen that the effects of geometric dispersion 
disappear for f/x = n or (1 — f)/x = m where n and 
m are integers. This corresponds to regions in which 
the scattered waves interfere destructively in the back- 
ward direction and the dispersion becomes the trivial 
one: KA = k a a + k^b. In these regions, the Bloch-wave 
is essentially a plane wave (i.e. it has only one dominant 
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FIG. 4: Dispersion surfaces for (Top, Left): TE and (Top, 
Right): TM Bloch waves with frequencies CU1...5. (Bottom, 
Left) : Dispersion surface for the leading TE wave components 
of waves CJ1...3. (Bottom, Right): The dispersion surfaces for 
TE and TM waves shown together. In the long wavelength 
limit the surfaces resemble those of a uniaxial birefringent 
crystal. As the wavelength is decreased to a length compa- 
rable to the periodicity of the structure, the surfaces increas- 
ingly distort. 



Fourier component). 

Figure [3| (Left) shows the electric field amplitude and 
(Right) the Fourier components of Bloch waves that cor- 
respond to the dots in Figure [3 As we approach the 
band edges the Bloch waves can be seen to differ strongly 
from plane waves. At the edge of a stop-band the ^2(^3) 
waves take the form of standing waves with the electric 
field concentrated in the material of higher (lower) refrac- 
tive index. The magnitude of K is the same at these 
two frequencies. When estimating the down-conversion 
amplitude it is useful to be able to switch between con- 
figuration and Fourier space depending on the degree of 
localization of the Bloch waves. For waves with non-zero 
fc|| both the interaction strength and the position of the 
bands will change. In particular the difference between 
TE and TM waves will become increasingly marked, as 
can be seen by inspection of the band structure of Fig- 
ure [21 Some intuition on how the magnitude of the Bloch 
momentum of propagating waves varies as a function of 
ku can be gained by examining Figure 21 where we show 
the dispersion surfaces for waves of frequencies CJ1...4. In 
particular it is apparent that in the long wavelength limit, 
the surfaces resemble those for a uniaxial birefringent 
crystal. The analogy can be formalized 13] by expanding 
the expression for K(kn, uo) fEa. 121)1 in the limit A ^> A to 
obtain expressions for the corresponding ordinary (TE) 



FIG. 5: Band diagram for the example Alo.4Gao. 6 As/Air 
structure with periodicity A = 187. 5nm. Points in the grey 
bands correspond to propagating TM (Left) and TE (Right) 
states. Points in the white bands correspond to frequencies 
that cannot propagate in the structure. The shaded area be- 
low 640nm represents the region in which Alo.4Gao.6As is ab- 
sorbing. 



and extraordinary (TM) refractive indices n and n e 



A 

a 

A" 



bl_ 

An? 



(26) 
(27) 



For phase-matching between waves with different polar- 
ization in the long- wavelength limit, optimization of the 
birefringence can provide a method to choose the fill frac- 
tions of a structure. 

Having gained some intuition on how the properties 
of Bloch-waves depend on the photonic-crystal geometry 
and constituent materials we now turn to the problem of 
calculating the down- converted emission from real struc- 
tures. 



VI. THE PHASE-MATCHING PROBLEM IN 
REAL STRUCTURES 

The parameter space for the design of a structure that 
achieves a particular type of phase- matching is large. 
Here we give a procedure to make predictions for the 
emission. We suggest iterating the procedure, making 
changes to the structure guided by the intuition devel- 
oped in the previous section, as a means of optimization. 



8 



750nm 










0^2 


































n 9 
















)i 




■ wm 












3 
































































0.0 


03 O.C 


22 


—on 


T7 

□ o.c 


05 O.C 


01 



FIG. 6: (Left): Field strength and (Right): Fourier compo- 
nents of degenerate TE/TM Bloch waves propagating along 
the axis of symmetry of the example structure. 



We take as an example a structure composed of 30 al- 
ternating layers of Alo.4Gao.6As (123nm thick) and air 
(64.5nm thick) and consider down-conversion of 750nm 
photons to degenerate 1500mn photon pairs. The frac- 
tion of Aluminium was chosen to avoid absorption of the 
pump photons. The relative thickness of the layers was 
chosen to optimize the birefringence in the long wave- 
length limit. The periodicity was chosen to ensure both 
pump and down- converted photons would propagate, the 
pump in the second band, the down-converted photons 
in the first band. We proceed according to the following 
recipe: 

(i) Band structure: We first plot the band diagram 
(FigureJSJ) of the structure. An intuitive variable in terms 
of which to plot the band diagram is the angle of prop- 
agation of the plane- wave solution inside the material of 
lower refractive index which is given by: sin(0) = 

When | ^ | > 1 total internal reflection occurs at the in- 
terfaces between the layers and the Bloch waves do not 
propagate along the crystal, but rather in a direction per- 
pendicular to the crystal axis. In the case of the example 
structure this occurs when k\\ is equal to the magnitude 
of the free-space wave- vector ~. It is of course neces- 
sary to include natural absorption of the materials in 
the band diagram plot since it can restrict the accessi- 
ble bands. Alo.4Gao.6As absorbs at wavelengths below 
640nm [l^. The natural dispersion of the materials is 
accounted for by including the frequency dependence of 
n ayb (See Ref. [l| for AlGaAs) in Eq.EEl 

(ii) Fourier spectrum: Having chosen frequencies 
that propagate (or having varied the periodicity to en- 
sure they propagate), we look at the amplitude of the 
Fourier components of the Bloch- waves having k\\ =0. 
Figure |H1 shows the waves and Fourier components for our 
example structure. Both have a leading Fourier compo- 
nent. The phase- matched emission will in general not 
have fen =0. However if the combination a;, fen =0 does 
not lie close to a band edge, then the fen = spectrum is 
a good indication of what the k\ \ 7^ spectrum will look 
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FIG. 7: Diagram to determine whether down-converted light 
will phase-match. The solid (dashed) blue lines represent the 
dispersion surfaces for pump (A = 750nm) TE(TM) photons. 
The solid (dashed) red lines represent the dispersion surface 
for down-converted (A = 1500nm) TE(TM) photons. The 
pump dispersion surface is centered at the origin, and the 
blue arrow represents the pump Bloch- vector. Two down- 
converted photon dispersion surfaces are drawn, centered on 
the extremities of the blue arrow. Intersections between the 
down- converted photon surfaces represent Bloch vectors that 
phase- match according to Ea ll7l 



like. 

(hi) Phase-matching of leading terms: We then 
look to satisfy the phase-matching equation K p = Ki + 
K2 + G starting with the leading G = n^z's from 
the Fourier analysis above. This is done by drawing a 
dispersion diagram such as the one shown in Figure [3 
Intersections between the displaced down-converted pho- 
ton dispersion surfaces represent Bloch vectors that 
phase-match according to Equation IT7I 

The full set of solutions corresponds to the intersection 
of the full two-dimensional dispersion surfaces. However 
the existence of an intersection in the diagrams discussed 
above is a necessary and sufficient condition for the two 
surfaces to have an intersection. 

(iv) Plotting the emission: Having determined that 
the intersections exist, we look at the detailed direction of 
the emission. To do this analytically we make the simpli- 
fying assumption that the crystal transverse dimensions 
are infinite (L Xl L y ^> \/n a ,b)- This is both justified in 
most experimental situations and makes it possible to 
produce a manageable analytical prediction of the emis- 
sion. The transverse part of the phase-matching function 
then becomes a product of two delta functions that en- 
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sure: 



^x 



— U 2 ) 

^x 

Up) _ 1.(2) 



Substituting these relations into the longitudinal part of 
the phase-matching function (Eq. ^|with the ^-function 
replaced by the appropriate sine function) to eliminate 

either k^} y or fc£^, we can then plot the phase- matching 
function for each emission. The elimination procedure is 
equivalent to calculating a partial trace over the emis- 
sions of the twin photon. It should be noted that the 
partial trace calculation is much simpler than in the case 
of birefringent non-linear crystals such as BBO since the 
optic axis is by default aligned with the crystal axis. Fig- 
ure [SI shows such a plot for three types of processes: 
Type I where a pump photon with TM polarization 
down-converts to photons 1 and 2 with TE polarization: 
p(TM) l(TE) + 2(TE), Type II: p(TM) l(TE) + 
2(TM) and Type III: p(TM) l(TM) + 2(TM). The 
width of the emission rings corresponds to the allowed 
longitudinal momentum mismatch and is inversely pro- 
portional to the crystal length. As can be seen in Fig- 
ure |H1 it is an important parameter for the crystal size 
being considered here. For a discussion of the relation 
between the crystal thickness and the yield of entangled 
photon pairs in BBO, see fl7j . 

(v) Efficiency of the process: From the plots of the 
detailed emission, the efficiency for emission in directions 
of interest can be estimated. The values of k\\ for pho- 
tons 1 and 2 are read off the plot. The amplitude for the 
phase-matched process can then be calculated by per- 
forming the sum over the Fourier coefficients (Eq. ITB|) . 
or by numerical integration of the field amplitudes in the 
non-linear layers fEa lTS|) . In the case of the Fourier am- 
plitude calculation it will be necessary to use the coef- 
ficients x(G) from the Fourier expansion of x(r) which 



C (a) (r) = -E 1 



Xbe 



sin — n) 

v A 7 



+ x a e sin {-^ n ) 



The n = term is simply the weighted average of the x s: 

A 

For photons emitted in the directions corresponding to 
the intersection of the rings in Figure |BI (Middle) or Fig- 
ure there is one dominant combination of Gs: G x = 0, 
G p = jf-z and Gi = G2 = 0. For this combination the 
Fourier amplitudes are 0.66, 0.90, 0.99 and 0.98 respec- 
tively. The combined effect of the Fourier amplitudes 
leads to a factor 0.58. 

In addition, the tensor nature of the x^ interaction 
needs to be included. This is done by taking the con- 
traction of the unit polarization vectors of the relevant 
Fourier components of E p ,Ei,E2 with the x^ tensor: 
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FIG. 8: Three types of phase- matching in the 
Al0.4Ga0.eAs/Air structure. (Top) Type-I (TM^TE+TE), 
(Middle) Type II (TM^TE+TM), (Bottom) Type III 
(TM^TM+TM). The density plots in k\\ space represent the 
emission strength calculated using only the phase-matching 
function and do not include contributions from the tensor 
nature of the x^ non-linearity. The blue circle represents 
the values of k\\ that totally internally reflect between the 
layers. The Type III phase-matching is weak when compared 
to types I and II for this set of parameters and so was 
multiplied by a factor 10 to make it visible. 



X^ tensor of AlGaAs has 43m point group symmetry 
and has only three non-zero coefficients all having mag- 
nitude 200pm/V 5]. For the conventional (100) surface 
orientation, the crystalline axis coincides with the direc- 
tion of normal incidence that we have defined to be the 
z direction. For photons emitted in the directions corre- 
sponding to the intersection of the rings in Figure [5| (Mid- 
dle) or Figure this leads to a value of 0.53 x 200pm/V. 

We are thus in a position to compare the efficiency of 
the process to that in BBO. The overall down-conversion 
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VII. CONCLUSIONS 



-0.1 0.1 0.2 0.3 0.4 0.5 

kx {n I A) 

FIG. 9: Polarization-entangled photon pairs can be extracted 
from the intersection of the TE and TM rings of the Type II 
emission according to a well-known scheme [j, El • 



efficiency for our example structure and choice of pump 
parameters is (0.53 x 0.58 x 200) 2 ~ 780 times that in a 
BBO crystal of similar size, where we have assumed that 
there is no reduction of the 2.2pm/V value of the BBO 
X^ 2 ) coming from the tensor nature of the interaction. 

It should be noted that there will also be down- 
conversion into directions that phase-match but do not 
involve leading Fourier amplitudes. These can be treated 
in exactly the same way, but will have much lower am- 
plitude. 

(vi) Extraction of polarization-entangled photon 
pairs: To extract entangled photon pairs, we proceed 
as in the scheme for the extraction of entangled pho- 
tons from non-linear crystals 3], by collecting photon 
pairs from the intersection between the TE emission and 
the TM emission (see figure 0J). At each intersection the 
photons will be polarized TE(H) or TM(V) with equal 
amplitude. However since the photons are emitted in 
TE/TM pairs, which ever the polarization of a photon at 
one intersection, we know the twin photon at the other 
intersection will have opposite polarization. Adding the 
amplitudes we obtain a maximally entangled polarization 
state: 

-L(\H)\V)+e ir \V)\H)) (28) 

where T is a phase that can be easily tuned experimen- 
tally, for example by placing a birefringent element in 
one of the paths. 

The existence of a fixed phase relation between the 
two terms depends on the genuine impossibility of de- 
ducing the polarization of a photon collected at one of 
the intersections from any of its other properties. Such a 
coupling between degrees of freedom reduces the coher- 
ence between the terms. There are many techniques to 
recover the coherence, for example the use of compensat- 
ing crystals Q or of a polarizing beam-splitter pj| UJ1 . 



We have presented a fully quantum-mechanical treat- 
ment of the down-conversion process in non-linear pho- 
tonic crystals, showing how to calculate the emission an- 
alytically in one dimensional structures. We have applied 
the theory to a realistic one-dimensional structure that 
consists of alternating layers of Alo.4Gao.6As and Air and 
demonstrated that entangled photons with a wavelength 
of 1500nm can be extracted from the down-conversion 
emission that results from the decay of pump photons 
having a wavelength of 750nm. 

We now suggest some possible extensions of the present 
work. The results for the mediated emission of sec- 
tion IVII are valid for structures of one, two and three 
dimensions. An obvious extension is thus the explicit 
analysis of Bloch-wave phase-matching in two and three- 
dimensional structures. An advantage of the use of two- 
dimensional structures is that two-dimensional structures 
with a large index contrast can be more easily fabricated 
than one dimensional structures with a large index con- 
trast. Another interesting extension is the analysis of 
down-conversion in the case where the pump photons ar- 
rive in short pulses or wave-packets (see Ref. 20] for a 
discussion of the problem in BBO). This brings up the 
topic of group velocity dispersion in photonic crystals 
which can differ significantly from its counterpart in nat- 
ural crystals. It seems likely that the possibility of tuning 
this form of dispersion will provide interesting possibili- 
ties. Another interesting avenue of research might be to 
investigate the potential advantages of down-converting 
to frequencies near a photonic band edge, where the den- 
sity of states is considerably larger than that at the centre 
of a band, thus further improving the efficiency. Finally 
it might be interesting to consider down- conversion of 
frequencies close to the non-linear material absorption 
band-gap, where the non-linearity can be significantly 
enhanced. 

In conclusion, the high degree of control afforded 
by the freedom of choosing the crystal geometry com- 
bined with the high non-linearity of semiconductor ma- 
terials for which many advanced fabrication techniques 
have been developed make non-linear photonic crystals a 
promising source of entangled photon pairs. 



APPENDIX A: DERIVATION OF THE 
BLOCH-WAVE AMPLITUDES 

In this Appendix we derive an expression for the 
Fourier coefficients of the periodic part of the electric 
field in one-dimensional Bragg structures. Throughout 
we will drop the TE(M) labels since the derivation is 
identical for the two cases. We begin by combining Equa- 
tions [221 and f° r the eigen-amplitudes of the forward 
and backward moving components of the electric field in 
each layer, with the expression for the electric field in the 
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nth layer (Eons. H$l and to obtain: 

E a = e iM-K.r) ^ (a +.a )e -ifeJ(z-nA) e iK 2 (z-„A) + 



\ ik^(z — nA) iK z (z— nA) 



(Al) 



Kr) ( (6+4)e"^ (2 "" A) e iK ^ (z - nA) + 



6 \ p ifeJ(2-nA) JK 2 (2-nA) 



.(A2) 



Comparing these relations with Eq. EI then yields the fol- 
lowing relation: 



n layers 

+ < n (z)4+<„(z)£ b _ (A3) 

Where E^_(z) represent the electric field propagating 
in the forward/backward direction in material a(b). In 
the nth layer they are given by: 

E ±,n( z ) = 4e^^- nA V K ^- nA ) (A4) 

4 >n W = b^e Tik ^ z - nA) e iK ^ z - nA \ (A5) 

To find the amplitudes of the of the Fourier coefficients 
sk,a we evaluate the Fourier transform of both sides of 
Eq. IA3I We define the Fourier transform / of a function 
/ as: 

1 



f(q) 



dzf(z)e~ iqz 



The Fourier transform of the left hand side of Ea. lA3l is: 

2tt> 



2tt\ 



(A6) 



To evaluate the Fourier transform of the right hand side 
of Ea. lA3l it is convenient to first re-express the sum over 
layers in terms of the following: the top hat function 



T XuX2 (x) - | Q ' x \ 
with X2 > xi, the comb function II(A,2iV-i-i) ( z ) : 



< X < X2 
< X\,X > X2 ' 



N 



n 



(A,2iV+l) 



(z)= E S(z-n\), 



-N 



E±(z) and E±(z) are well defined functions at all points 
in the the crystal. They are equal to outside layers 
a and b respectively and equal to E± n (z) (Eg. IA4|) and 
E± n (z) fEa. lA5|) in the nth layer. The Fourier transform 
of E±(z) and E±(z) can be simplified using the standard 



relation: 



to obtain: 



E± 
El 



f*g = f~g (Aio) 

(^±,n=0 T -a,o) n (A,2AT+l) 
(^l n ^A,-a)n ( A,2iV + l). (All) 



If we evaluate the Fourier transform of II(A,2iV+i) : 

1 sin [q4(2N+l)]) 
[n (A , 2jV+1) ]fe) - -j= — ^ , (A12) 



2tt 



and take the limit of N — > oo, we obtain: 



N 



lim II( Ai 2jv+i) = v27rll 



(A13) 



comparison in this limit with the left hand side of Eq. IA6I 
then leads to the following expression for Sk,a : 

*k,a(»x) = 

E<^fl afi (n^-)e a + + Ef^C afl (n^)e a _ 



Evaluating E% T- afi (q z ) and E b ±fi T-^- a (q z ): 
E<Qfl afi (q z ) = 



u (K zT k a z -q z )% 

E b ± ^lC A ^ a (q z ) = 

: sm[(K z Tk b z -q z )l] 



(A15) 



V 7 ^ U (tf* T - q z ) 



and the convolution of functions / and g: f * g: 

/oo 
dzf(z)g(x- z). 
-co 

In terms of these functions we have: 

Eiayers KA^l + ^-,n(^- + + 

£+(z)£+ + ££(«)e^ + £+(«)e+ + E b _(z)£ b _. 



completes our derivation. Combining Eq. IA15I and 
Eq. rUTil with Eq. EH 

we obtain the result given in 

section Ivl 
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